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ABSTRACT 

Polarimetric studies of pulsar radio emission traditionally concentrate on how the 
Stokes vector (/, Q, V) varies with pulse longitude, with special emphasis on the 
polarization angle (PA) swing of the linearly polarized component. The interpretation 
of the PA swing in terms of the rotating vector model is limited by the assumption of 
an axisymmetric magnetic field and the degeneracy of the output with respect to the 
orientation and magnetic geometry of the pulsar; different combinations of the latter 
two properties can produce similar PA swings. This paper introduces Stokes phase 
portraits as a supplementary diagnostic tool with which the orientation and magnetic 
geometry can be inferred more accurately. The Stokes phase portraits feature unique 
patterns in the /-Q, I-U ^ and Q-U planes, whose shapes depend sensitively on the 
magnetic geometry, inclination angle, beam and polarization patterns, and emission 
altitude. We construct look-up tables of Stokes phase portraits and PA swings for pure 
and current-modified dipole fields, filled core and hollow cone beams, and two empirical 
linear polarization models, L/I = cosOq and L/I = smOi^, where 6*0 is the colatitude 
of the emission point. We compare our look-up tables to the measured phase portraits 
of 24 pulsars in the European Pulsar Network online database. We find evidence in 
60% of the objects that the radio emission region may depart significantly from low 
altitudes, even when the PA swing is S-shaped and/or the pulse- width-period relation 
is well satisfied. On the other hand, the data are explained adequately if the emission 
altitude exceeds ^ 10% of the light cylinder radius. We conclude that Stokes phase 
portraits should be analysed concurrently with the PA swing and pulse profiles in 
future when interpreting radio pulsar polarization data. 

Key words: magnetic fields — polarization — pulsars: general 



1 INTRODUCTION 

Polarimetric studies of pulsar radio emission probe the ge- 
ometry of the pulsar magnetosphere. Although the phe- 
nomenology of the emission is often complex, there are some 
general trends in the observed pulse shapes and polarization 
angle (PA) swings that have been explained qualitatively 
with simple magnetospheric models. In many objects, the 
data are approximately consistent with low-altitude emission 
from open dipolar field lines centred on the magnetic axis 
(e.g. lRankinlll983 : iLvne fc ManchesteJlToSSal '). although in 
some objects, like the Crab, the emission is thought to 
emerge as a fan be am from higher altitudes in the outer 
magn etosphere (e.g. ICheng fc Rudermanlll977l : iKaspi et al.l 
120001 ). 

Three lines of evidence favour an approximately dipolar 
magnetic field in the radio emission region. First, pulses nar- 
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row with increasing frequency. According to the radius-to- 
frequency mapping, the frequency of the radiation is related 
to the local particle density, whic h decreases with altitude, 
as th e dipole field lines diverge (|Ruderman fc Sutherland! 
1 19751 : ICorded Il978l 1. The width of the pulse profile is 
commonly used to infer the emission height, with high- 
frequency p u lses o rig inating closer to the surfa ce (e.g . 
Manchesteil 1 19961: Ivon Hoensbroech fc XilourisI 1 19971: 



Karastergiou fc JohnstonI 20071 : Weltevrede fc JohnstonI 
2008h . Second, in many objects, the pulse width W 



is inversely related to the spin period P of the pulsar 
according to oc P"^/^ (|Rankinl [l993 ). This relation 
can be explained geometrically assuming an emission 
cone tangent to the last op en field line of a dipole (e.g. 
iGangadhara fc Guptal l200lh . Under this assumption, the 
measured pulse widths can also be related to the cone 
half-opening a ngle and the ang le between the magnetic and 
rotation axes (|Gil et al.l Il984 ). Third, the plane of linear 
polarization traces out a characteristic S-shaped swing over 
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one pulse period. The swing has been modeUed successfully 
by the rotating vector model assuming a dipole field (or, 
more generally, an axisymmetric field). The model relates 
the gradient of the S-shape at the infiection point to the 
inclination angle of the obse rver and the angle between the 
magnetic and rotation a xes dR adhakrishnan fc Cook3ll969l : 
iLvne fc Manchestei|[T98 8b: Hib schman fc Aronsll2o'oif ). The 
phase lag between the centroid of the pulse profile and the 
infiectio n point has also been use d to estimate the emission 
altitude (|Blaskiewicz et al]ll99lh . 

Unfortunately, experience shows that there are limita- 
tions in relying only on the pulse profile and PA swing. For 
example, the three trends mentioned above fail to fix the ge- 
ometry uniquely. The PA swing looks remarkably similar for 
many different orientations of the observer and the magnetic 
axis. Moreover, observations show that many pulsars do not 
display the clean, S-shaped swing expected from the rotat- 
ing vector model. Instead, the S-shape is often distorted, 
suggesting nondipolar con f igurations (iKarastergiou et al.l 
I2OO5I : iJohnston et ai]|2008al : iHan et al.|[2009l ). I n particular 
many millisecond pulsars have fiat PA profiles (|Stairs et al.l 
1 19991 : lOrd et all [2004.'). highly d istorted profiles (e.g. PSR 
J0437-4715: iNavarro et"al]ll997l ') or seemingl y random PAs 
(e.g. giant pulses from PSR J1824-2452A; iKnight et al] 
I2OO6I ). Multi- frequency observations show th at the PA swing 
varies with frequency for many pulsars (|johnston et al.l 
l2008bl ). indicating that the magnetic geometry changes ap- 
preciably with altitude. Several mechanisms can distort a 
dipole field, e.g. a current flowing along the field lines 
l|Hibschman fc Aronsll200ll : lDvkdl2008l'). or rotational sweep- 
back near the light cylinder ([Hibschman fc Arond I2OOII : 
iDvks fc Hardinell2004l : lDvksll200^ 

One diagnostic tool which has been used sparingly but 
holds considerable promise is the Stokes phase portrait, that 
is, the pattern traced out by the four Stokes parameters /, 
Q, U and V when plotted against each other over one pulse 
period. PA swings describe how the quantity U/Q varies 
with time, but less effort has been devoted to studying how 
Q and U (say) depend on each other. We show here that 
Stokes phase portraits harbour a great deal of extra infor- 
mation regarding the magnetospheric geometry. The analy- 
sis of polarization data on the Q-U plane is common prac- 
tice when calculating rotation measure transf er functions 
for polarized emission from radio galajdes (e.g. Burnlll966l : 
iBrentiens fc de Bruvn|[2005l : iHaverkorn et al.ll2006l). and in 
optical and UV pulsar data analysis (e.g. Smith et al.ll 19881 : 
iGraham-Smith et al.lll996l : ISlowikowska et al.ll2009l ). 

In this paper, we calculate theoretically how the Stokes 
parameters vary with pulse longitude for several magnetic 
field configurations and beam patterns. We apply the the- 
oretical results to several radio pulsars to determine the 
observer's inclination angle, i, the angle between the mag- 
netic and rotation axes, a, the emission altitude, and the 
approximate geometry of the magnetosphere (e.g. the ra- 
tio of poloidal to toroidal field). We call this reconstruction 
technique Stokes tomography. In Section [2] we define the 
model and the algorithm used to generate the Stokes phase 
portraits. Look-up tables of phase portraits, PA swings, and 
pulse profiles for a pure dipole field at low emission altitudes 
are presented in Section |3] The look-up tables can be com- 
pared directly against observational data for any pulsar to 
determine the emission geometry. In Section [l] we present 



look-up tables for a pure dipole at 10% of the light cylinder 
radius, where relativistic aberration effects are important. 
It is shown that many pulsars that are categorised as pure 
dipoles on the strength of their PA swing possess Stokes 
phase portraits that are inconsistent with a pure dipole at 
low emission altitudes. In Section O look-up tables are pre- 
sented for a current-modified dipole field with a radially in- 
creasing to roidal component, as in s tandard magnetospheric 
models (jHibschman fc AronsI I2OO1I ') . Future applications of 
Stokes tomography are canvassed briefiy in Section [6] Com- 
panion papers will generalize the approach to circularly po- 
larized pulsar radio emission and realistic (e.g. force-free) 
magnetic fields. 



2 STOKES TOMOGRAPHY 
2.1 Radiation field 

In the magnetosphere, radiation from highly relativistic par- 
ticles tied to magnetic field lines is narrowly beamed along 
the particle velocity vector. The observed emission point 
xo(t) at any time t is therefore located where the mag- 
netic vector B[Xp(t),(t)] points along the observer's line- 
of-sight vector w. The unit tangent vector to the magnetic 
field at xo{t) is defined as t = B[xo(t), t]/|B[xQ(t), For 
non-relativistic particles, to find the emission point at time 
t, we must therefore solve the equation t = w. In general, 
however, when considering emission altitudes of r > O.lrLC, 
where tlc = c/fi is the light cylinder radius, we must ac- 
count for relativistic aberration. Aberration shifts the emis- 
sion point and electric field vector by different amounts of 
order r/r ^c- T he emission point xo(t) at time t satisfies the 
equation (|Blask icwicz ct al. 1991) 



t + flxxo/c 
\t + fl X xo/c| 



(1) 



where Jl is the angular velocity vector. Note that B varies 
with t intrinsically, as the star rotates in the observer's 
frame, not just through xo(t). The normal to the field at 
xo is deflned as n = k./|k|, where k — {t ■ V)t is evaluated 
at xq. The binormal is defined as b = t x n. 

The amplitude and direction of the complex electric 
field vector E at the emission point are determined by the 
radiation physics, which is not understood in detail. A pop- 
ular assumption is that E at time t points along the instan- 
taneous acceleration vector of the par ticle, a. To first order 
in r/r-Lc, the acceleration is given by l|Dvksll2008D 

dt 

a = c— + c(v ■ V)t + cr2xt + nx(nx Xo), (2) 
dt 

in terms of the instantaneous velocity vector v — ct + f2 xxq, 
assuming no precession {dft/dt = 0). From left to right, the 
four terms on the right-hand side of ((2} describe: (i) the 
change in field aligned velocity ct as the magnetic dipole 
moment (and hence t) rotates infinitesimally in the iner- 
tial frame; (ii) the change in field-aligned velocity ct as the 
charge moves infinitesimally along v = ct-|-i7 x xq and hence 
t changes; (iii) the change in corotation velocity 17 x xo as 
the charge moves infinitesimally along the field t; and (iv) 
the change in corotation velocity f2 x xo as the charge coro- 
tates infinitesimally with the magneto sphere along x xo 
(|Dvksll2008l : iHibschman fc Aronsll200lh . 
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The associated Stokes parameters (7, Q, U, V) describ- 
ing the polarization state are defined as 



/ = \E^\^ + \Eyf (3) 

Q = \E^\^^\Eyf (4) 

U = 2Re{E,E;) (5) 

V = 2Im{E^E*y), (6) 



where / is the total intensity, L = (Q^ + C/^)^/^ is the hn- 
early polarized component, and V is the circularly polarized 
component. The x- and y- components are measured with 
respect to an orthonormal basis (x, y) which is fixed in the 
plane of the sky. In this paper, we choose x = flp/lilpl and 
y = X X w, where f2p = f2 — (f2 ■ w)w is the projection of 
the angular velocity CI onto the plane of the sky. Then the 
polarization angle, tp, measured counter-clockwise between 
X and the linearly polarized part of E is given by 

V.= -tan -. (7) 
2.1.1 Orientation of i and y 

In the observational data, the direction of fl is usually 
not known. The polarization basis in which the observed 
Stokes parameters are expressed may therefore differ from 
the canonical basis. This affects the Stokes phase portraits 
in the following way. 

Consider two polarization bases, (x, y) and (x', y'), 
where (x', y') is rotated clockwise by an angle /3 with respect 
to (x, y). In the canonical (unprimed) basis, the electric field 



vector transforms as 

£^ = cos l3E'^ + sin l3E'y, (8) 

Ey = - sin 13 E'^ + cos pEy, (9) 

and the Stokes parameters transform as 

I ^ I' (10) 

Q = Q' cos 2/3 + (7' sin 2/3 (11) 

U = -Q' sin 2/3 (7' cos 2/3. (12) 



Hence, in the Q-U phase portrait, the pattern rotates 
counter-clockwise through an angle 2j3 with respect to the 
Q'-U' phase portrait, but its shape remains invariant. In 
contrast, in the I-Q and I-U phase portraits, both the shape 
and orientation of the pattern change. When fitting our 
model to data, in cases where /3 is unknown, we first work 
to match the shape of the pattern in the Q-U data, without 
worrying about the orientation. Then, armed with a good 
fit to the Q-U shape, we rotate the basis to match the orien- 
tation as well and deduce /3. Finally, knowing P, we match 
the patterns in the I-Q and I-U data. Note that the PAs in 
the two bases are related hy ip = ip' ~ P , as expected. 

2.2 Numerical algorithm 

We now describe how to construct (J, Q, U, V) as functions 
of pulse longitude (and hence each other) for a prescribed 
magnetic geometry. 

Let us define two reference frames, as depicted in Fig- 
ure [T] the inertial frame in which the observer is at rest, 
with Cartesian axes (e^,, 6^,63), and the body frame of the 



rotating pulsar, with Cartesian axes (ei, 62, 63). The line-of- 
sight vector w is chosen to lie in the Sy-e^ plane, making an 
angle i with e^. The rotation and magnetic axes are chosen 
to lie along ez and 63 respectively. The orientation of the 
body axes with respect to the inertial ajces is computed as a 
function of time by solving Euler's equations of motion for a 
rigid body; precession can be included in general, although 
we do not examine it in this paper. 

We define a spherical polar grid (r, 6, <j)) in the body 
frame covering the region Xmin ^ r/r-^c ^ a;max, ^ 6* ^ tt, 
^ (/> ^ 271, with 64 X 256 x 64 grid cells, where the line 6^0 
lies along 63. Most of our analysis is done in the middle to 
outer magnetosphere, with a;inin = 0.001 and a::max = 0.38. 
Nonetheless, other choices are possible and easy to imple- 
ment. Given a magnetic field B(r, S,(/>), we compute t, n, 
and b at each grid point using second-order central finite 
differencing and ghost cells at the edges. 

In order to calculate (J, Q, U, V) versus t, we search the 
grid at a fixed altitude ro to find the emission point satisfy- 
ing (m (the point P in FigureQ. At P, (t 4- H x xo/c) ■ w is 
a maximum; note that the location of P chariges with time 
in both the body frame and the inertial frame|j To fine-tune 
the location, we interpolate t linearly between the four sur- 
rounding points and re-compute ([1} for each interpolated 
point. In general, there are two points xo(t) satisfying ((T} 
on opposite sides of the pulsar, but the observer only sees 
the one facing him. 

We stipulate the shape of the beam pattern I{r,9,(j)) 
and the degree of linear polarization L{r,6,4)), and assume 
that E is constant and lies along a, i.e. 

E = ia/7. (13) 

In Figure [T] I{r, 8, (j)) is represented by the shading on the 
surface of the sphere. It peaks along 63 (i.e. ^ = 0), i.e. 
Figure [1] depicts a centre-filled beam. Other beam patterns, 
e.g. a hollow cone, are also considered in Sections [3] and [S] 
Note that E is purely real in this paper because we assume 
V — 0, i.e. zero circular polarization. Examples of I and 
L are given in Sections 12.41 and 12.51 drawn either from a 
phenomenological theory or directly from data. 

2.3 Relativistic aberration 

Relativistic aberration acts to shift the pulse centroid rela- 
tive to the inflection point of the PA swing, such that they 
are separated in phase by » 4r/rLC radians to leading order 
(Bl askiewicz et al.l llQQll : iHibschman fc Aronj I2OOII : iDvksl 
I200C Hence, aberration introduces an altitude dependence 
in the case of a pure dipole, which is absent in the non- 
aberrated case. At large emission altitudes (> O.lr-Lc), the 
phase shift can significantly alter the shapes of the Stokes 
phase portraits. This is investigated further in Section U) 
Additionally, if the magnetic field sweeps back rotationally, 

^ In general, we must consider ±t in JTJ to capture emission 
from particles flowing out from both north and south poles. In 
this paper, however, we restrict out attention to emission from 
one pole only, to avoid cluttering the look-up tables in Sections 
[3] and [5] We therefore do not capture "interpulse" emission with 
our model. The 24 objects studied in Section 15.31 exhibit either 
zero or weak interpulse emission in existing data. 
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ei' 

62 

Figure 1. Orientation of the inertial axes ex,ey,ez (solid lines) 
and pulsar body axes ei , 62 , ea (dashed lines) at an instant where 
62 and 63 lie in the 6^-6^ plane. The line of sight vector w 
(dashed-dotted line) is defined to always lie in the Sy-Bz plane 
at an angle i to 62. 63 is tilted with respect to the rotation axis 
Bz by an angle a. An example of a dipolar field line (B) and its 
corresponding emission point xo(t) (labelled P) are drawn. At P, 
the tangent vector is t = w and the normal vector is n. The an- 
gle between t and 63 is 60 • The shading on the sphere represents 
the beam pattern I{r,9,<p); I is a maximum where the shading is 
brightest. 

or develops a helical, current-modified component, the emis- 
sion point at an instant in time varies as a function of al- 
titude. The effects of aberration combined with a current- 
modified field are investigated in Section [S] According to 
the radius-to-frequency mapping, observing emission from 
different altitudes is analogous to making observations at 
several different frequencies. 

Equation ((T]) ignores cross-field particle motions (e.g. 
curvature and gradient drifts) in the outer magnetosphere, 
where the cyclotron cooling time exceeds the fiow time, for 
the sake of simplicity. 

2.4 Beam pattern 

The two most common be am patterns in the literature are 
the core-and-cone model ('Rankin 'l983') and the patchy- 
beam model (Lyne & Manchester 1988a). In the core-and- 
cone model, single-peaked profiles are explained by emission 
filling the polar fiux tube, double-peaked profiles are ex- 
plained by emission from a hollow cone tangential to the 
polar fiux tube, and triple-peaked profiles are explained by 
a combination of the two. The patchy-beam model was pro- 
posed to explain asymmetric profiles. It suggests that emis- 
sion comes from randomly distributed locations within an 
emission beam. In this paper, for simplicity, we consider two 
beam shapes based on the Rankin model: (i) a filled core 



beam (single-peaked pulse profile) and (ii) a hollow cone 
(double-peaked pulse profile). Both beams are centred on 
e3 and are cylindrically symmetric about it. 

The colatitude of the emission point at time t, 6o{t), is 
obtained from 60 = cos~^[r^^xo(t) ■ e3] and can range over 
^ ^0 ^ TT- We model the polarized intensity I{t) = /[So(i)l 
as a gaussian beam, 

I{eo) = (27r<7^)-^/"exp [-(^o - pf/i^^^)] (W) 

where a defines the beam width, and p is the half-opening 
angle of the emission cone. In general, a and p vary from 
pulsar to pulsar. However, to generate the look-up tables in 
Sections[3HSI we use fixed, representative values: p = 0° and 
a — 10° for a filled core, and p = 25° and a = 10° for a 
hollow cone. For the sake of simplicity, we assume in what 
follows that only one pole shines (see footnote above); it is 
trivial to modify this assumption if necessary. 

Equation p4[) can only be used to model approxi- 
mately symmetric pulses. To fit asymmetric profiles, more 
complex models, such as a 'broken cone' or 'horseshoe', 
are required, in which the intensity varies with longitude 
(poit) as well as 9o{t). Such asymmetri c intensity maps are 
inspired by the patchy beam m odel (|Lvne fc Manchester! 
ll988al : lKramer fc Johnstonll2008l ). 



2.4-1 Tilted axis of symmetry 

Although the magnetic-pole model for pulsar radio emission 
is commonly adopted, it is possible that the beam pattern is 
centred on another axis that is slightly tilted with respect to 
the magnetic axis. Fan beams in the outer m agnetosphere 
arc an extreme instance of such a geometry fehen g et al.l 
|200_Q; _Wattcrs et al. 2009.), but less dramatic displacements 
can occur if there are local, quadrupolar surface fields in the 
vicinity of the magnetic poles. 

To illustrate the above effect, we consider an arbitrary 
example, where I and L are symmetric about the axis e3 in a 
reference frame (e'j, e2, e'^), where e3 is tilted with respect to 
e3 by 10° , and e'l is tilted with respect to ei by 90° . Figure[2] 
compares the original (top panel) and tilted (bottom panel) 
pulse profiles, PA swings, and Stokes phase portraits for 
(a,j) — (30°, 40°). The axis tilt causes the pulse phase to 
lead the PA swing by f» 0.47 radians. In the bottom panel, 
the hockey stick in the I-Q phase portrait broadens into a 
balloon, while the shape in the I-U plane narrows and tilts. 
The heart shape in the Q-U plane tilts and is no longer 
symmetric about U = 0. 

This effect mimics relativistic aberration, which results 
in a relative phase shift between the pulse profile and PA 
swing of 4r/rLC rad. The effects of relativistic aberration on 
the phase portraits are investigated in detail in Section|4l We 
note that neither Stokes tomography nor the rotating vector 
model can distinguish the two effects. For the sake of clar- 
ity, in Sections [SHS] of this paper, we therefore assume that 
the beam and polarization patterns are symmetric about e3, 
while noting that the Stokes tomography technique can be 
generalized easily (as in Figure [2]) to accomodate a tilted 
axis of symmetry. When modelling individual pulsars in de- 
tail, / and L can be fitted iteratively and empirically, and 
asymmetric beam patterns arise naturally. 
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Figure 2. Theoretical pulse profiles, PA swings, and Stokes phase 
portraits for a pure dipole and a beam pattern symmetric about 
63 (top four subpanels), and about (bottom four subpanels), 
an axis tilted by 10° with respect to 63 and 90° with respect to ei, 
as described in the text. Within each half of the figure, the four 
subpanels display (clockwise from top left): (a) ///max (bottom) 
and PA (top; in radians) as functions of pulse longitude, (b) I-Q, 
(c) Q-U , and (d) I-U. Relativistic aberration is not included in 
this figure. 



2.5 Linear polarization 

The uncertainty surrounding the microscopic emission 
physics in a pulsar makes it difScult to model its po- 
larization properties, although theoretical progres s has 
been made (e.g. iGil fc Snakowski, ,1990.; ,Mehos3 Il995l: 



McKinnon fc Stinebrine 
20031 : 



19981 ; Rankin fc RamachandranI 



Xu et al.H200oi r Much of this work attempts to ac- 



count for observe d orthogonal mode .jump s (i.e. 90° jumps 
in the PA swing ). iMcKinnon fc Stinebrin"3 l|l998l ) modelled 
the L/I distribution statistically, describing Q as a Gaussian 
random variable. One physical model assumes that emission 
comes from inverse Compton scattering, and show s that the 
fract ional linear polarization L/I varies with Oo JXu et al.l 
I2OOOI ). 

In the absence of a unique empirical or theoretical rule 
for how the fractional linear polarization L/I varies across 
the beam, we investigate two simple cases when construct- 
ing the look-up tables in Sections [3] ancO (i) L/I oc cos 60 
(generally, but not exclusively, suited to single-peaked pro- 
files), and (ii) L/I oc sin So (generally, but not exclusively, 
suited to double component profiles). The proportionality 
constants are chosen to be unity in Sections [SHS] for sim- 



plicity, but, in reality, approximat ely 80% of normal p ulsars 
have a pulse-averag ed L/I ^ 0.3 (|Xilouris et al.lll998l ). 

Although the above models are useful for illustrating 
the variety of Stokes phase portraits in our look-up tables, 
they are a rough approximation in any specific object. Model 
(i) gives L/I 1 for small Oq, higher than the observed 
fraction in most pulsars. Moreover, in many objects, there 
is a phase lag between / and L, which we ignore. We do not 
attempt to model orthogonal mode jumps in this paper. 



2.6 Summary of algorithm 



Figure 3(a) summarises, in the form of a fiowchart, the steps 



required to generate the PA swing, pulse profile, and Stokes 
phase portraits for a pulsar with a given magnetic configu- 
ration (as described in Sections l2.1H2.5l above). Figure 3(b) 



summarises the steps required to match the theoretical PA 
swing, pulse profile, and Stokes phase portraits from the 
look-up tables to observational data and thereby infer the 
pulsar's orientation and magnetic geometry. For a nondipo- 
lar field, the results depend on altitude, so the process must 
be repeated using several values of r until a good match is 
achieved. We note that in principle, if /3 is unknown, one can 
also repeat the process using several values of p. However, 
throughout this paper, we find that once an approximate 
estimate for P is found, it is more useful to explore a and i 
in detail, as varying /? only alters the orientation of the Q-U 
pattern, and not its shape. 



3 DIPOLE FIELD AT LOW EMISSION 
ALTITUDES 

In this section, we calculate Stokes phase portraits, PA 
swings, and pulse profiles for a pulsar with a purely poloidal, 
dipolar magnetic field at low emission altitudes, where the 
effects of relativistic aberration are negligible. We generate 
a two-dimensional look-up table of phase portraits for ori- 
entations 10° ^ i < 90° and 10° < a 90°. For each beam 
pattern (core emission in Section [3. II conal emission in Sec- 
tion and polarization pattern (L/I = cos So, sin So), we 
present three phase portraits (l-Q, 1-17, Q-U) for each pair 
of angles (a,i). Finally, in Section [4.31 we analyse a selec- 
tion of full Stokes data obtaine d from the European Pulsar 
Networ k (EPNfl online archive (|Lorimer et al.lll998l ') for ob- 
jects in iGould fc Lv"n3 (|l998h . We show that a dipole field 
(and hence low-altitude emission) is inconsistent with the 
observed phase portraits of many superficially dipole-like 
pulsars, e.g. objects which obey the pulse-width-period re- 
lation or exhibit clean, S-shaped PA swings. 



3.1 Filled core emission 

3.1.1 L = /cosSo 

Examples of the / (solid curve) and L (dashed curves) pro- 
files for (a, i) = (70°, 20°) are shown in Figure |4l / is 
normalised by its peak value. The pulse is single-peaked 
and maximally polarized at the peak. The dotted curve 

^ Available at: http:/ /www.mpifr-bonn.mpg.de/pulsar/data/| 
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Construct observed phase portraits 
I'-Q', I'-U', Q'-U' from Stokes data 



Check look-up tables for closest matching 
shape in Q'-U' 

I 



Generate t, n, b in body frame (61,62,63) 



Rotate Q'-U' to match orientation correctly; 
measure /3 



Solve Euler equations for orientation 
of (61,62,63) in inertial frame (63,, 6^,62) 



Solve (1) to find emission point xo(t) 



Generate electric field vector at 
emission point, E = L/Ia. 

I 

Generate Stokes parameters using 
polarised components of E 
in canonical basis 



Generate PA swing from Q, U 
and phase portraits I-Q, I-U , Q-U 
in canonical basis 



(a) 



Use measured /3 to transform 
(/', Q' , U') data to the canonical basis (/, Q, U) 
using (12) and (13) 

i 

Generate Stokes phase portraits for a finer 
grid centred on closest matching (a, i) 

\ 

Use data to construct trial 
beam pattern I{B, (j)) for best-fitting (a, i) 



Use data to construct trial 
linear polarisation profile L{d, <p) 
for best-fitting {a, i) 



Compute Stokes phase portraits and 
PA swing using / and L 



(b) 



Figure 3. (a) Recipe to generate look-up tables of Stokes phase portraits, PA swings, and pulse profiles for a given orientation and 
magnetic geometry, (b) Iterative recipe to extract orientation and magnetic geometry from observational data. 



shows how the magnetic colatitude 8o{t) of the emission 
point varies across one pulse period. We note the follow- 
ing trends, (i) The pulse profile narrows with increasing i 
and Q as the emission point moves further from the mag- 
netic pole. The full-width half-maximum (FWHM) decreases 
from ~ 0.5 phase units at (a, z) = (10°, 10°) to ~ 0.15 at 
(a,j) = (90°, 90°). (ii) For a = i, where 60 = 0° at the 
pulse peak, we have L{t) — I{t). L{t) decreases as \a — i\ 
increases, e.g. we find L{t) « 0.6/(t) at (a, i) — (10°, 90°). 
The maximum value of L/I can be calculated in terms of a 
and i by solving equation (|15p below for cos6'o. 

In FigureO we plot the path traced by xo(i) in the body 
frame across one pulse period for three cases: (A) (a, i) — 
(20°, 10°), (B) (20°, 20°), and (C) (20°, 30°). For a dipole 
field, the emission colatitude 9o is related to a and i as 
follows: 



cos(j ■ 



3 cos^ Oo - 1 
(l + 3cos2 6lo)i/2- 



(15) 



The path changes from an undulation at a < i (curve C 
in Figure [SJ to an ellipse at a > i (curve A) . The differ- 
ence in geometry is largely responsible for the difference in 
polarization properties observed for a < i and a > i. 

The Stokes phase portraits for a filled core with L/ 1 = 
cos 6(i are drawn in Figures [GHS] In each figure, the panels 
are organised in order of increasing i (left-right) and a (top- 
bottom) in landscape mode. The panels range over 10° ^ 
i ^ 90° and 10° ^ a ^ 90° and are separated by intervals of 
10° . In Figures |S] and [T] / is plotted on the horizontal axis 
(normalised by its peak value) , and Q and U are plotted on 
the vertical axis. In Figure [S] Q is plotted on the horizontal 
axis and U is plotted on the vertical axis. The same layout 
is used in all the look-up tables presented in this paper. 

In the I-Q plane (Figure |6]), we note the following be- 
haviour, (i) For a ^ 30° and i ^ 50°, the pulse traces 
out an approximately straight line segment which evolves 
into a slender hockey-stick shape as ct increases and i de- 
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Pulse longitude, I 

Figure 4. Dipole field at a low emission altitude r <g tlc- Exam- 
ple of a pulse profile for a filled core beam with linear polarization 
L = /cos 6*0 and (a,i) = (70°, 20°). Solid, dashed and dotted 
curves represent the total polarized intensity /, degree of linear 
polarization L, and emission point colatitude Bq. Pulse longitude 
I is measured in units of degrees. 




(Ii (rod) 

Figure 5. Examples of the path traced by the emission point 
xo (t) in the body frame (curves) , overplotted on the intensity map 
for a filled core beam (greyscale). The beam is brightest at 9 = 
(corresponding to the es axis). Curves A: (a, i) = (20°, 10°), B; 
{a,i) = (20°, 20°), and C: (a,i) = (20°, 30°). 

creases. In the body frame, for a > i (below the diago- 
nal in Figure [6]), xo(i) traces out an ellipse in the quadrant 
0^9^ 7^/2, ^ TT (see Figure [5}. As \a — i\ increases, 
the ellipse moves further from 8 = 0° (and hence away from 
the beam peak). This reduces the curvature of the hockey 
stick. For a < i (above the diagonal in Figure|6]), xo(t) traces 
out a sinusoidal curve across ^5 (/> ^ 2tt. Again, the curve 
moves away from 6 = 0° as \a—i\ increases, reducing the cur- 
vature of the hockey stick. The left horn of the hockey stick 
shortens as xo(i) moves away from 6 = 0°. (ii) For a = i, the 
pattern is a straight line with dQ/dl < 0. In the body frame, 
xo(t) traces out a circle that crosses ^ = 0° symmetrically, 
(iii) The shapes are slightly larger in the panels adjacent to 
the a = i diagonal (i.e. -0.2 < Q < 1.2 for \a - i\ < 10°) 
and decrease to e.g. < Q < 0.6 at {a,i) = (10°, 90°) at 
the corner of Figure [B] 

In the I-U plane (Figure [7|), we note the following be- 
haviour, (i) The pattern is symmetric about U = and 



traces a sideways balloon shape for most angles, with the 
pointy tip at {U,I) = (0,0). (ii) For a s; 50° and i ^ 30°, 
the pattern twists into a figure-eight. This twist occurs when 
the path traced by xo(t) in the body frame elongates from 
a circle to an oval, (iii) For a = i, the phase portrait is 
roughly oval, with major and minor axes along the /- and 
?7-axes respectively. It narrows with increasing a = i. (iv) 
The patterns are broadest adjacent to the a = i diagonal 
(e.g. —1 < U < 1 for |a — i| < 10°) and narrow by up to 
90% as \a — i\ increases to the corners of Figure [T] 

In the Q-U plane (Figure (Sjl, we see a mix of ovals, bal- 
loons and heart shapes with the following properties, (i) For 
a < i, above the diagonal, the shapes broaden with increas- 
ing a, while the reverse is true for a > i, below the diagonal. 

(ii) On the diagonal, the ovals narrow with increasing a = i. 

(iii) In some cases, e.g. {a,i) = (50°, 40°), a cusp forms at 
the start and end points of the pulse, where ([/, Q) = (0, 0). 
This cusp twists into a secondary oval for 10° ^ a ^ 50° 
and i ^ 20°. As in the I-U plane, this secondary loop forms 
when the pattern traced by xo(f) elongates from a circle to 
an oval. The size of the loop is proportional to the ellipticity 
of the oval, (iv) As in Figure [T] the patterns in the phase 
portrait are broadest for |q — i| < 10° directly adjacent to 
the a = i diagonal, and narrow as \a — i\ increases. 

In Figure [9] we plot the PA swings corresponding to 
each panel in Figures [GHS] We plot only the parts of the 
swing that are illuminated by the pulse, i.e. when \L\ J5 
10^^. The threshold is arbitrary and should be adjusted 
for instrumental resolution when modelling actual data. We 
note the following behaviour, (i) The PA swings are iden- 
tical to those predicted by the rotating vector model, as 
expected for a dipole magnetic field, with the S-shape steep- 
ening as a and i increase, (ii) Above the a = i diago- 
nal, we find d{PA)/dl > 0, where / is the pulse longitude, 
whereas the opposite is true below the diagonal. This is 
consistent with the rotating vector model's prediction that 
d{PA)/dl = sina/sin(i — a) at the inflection point, (iii) 
Along and below the a = i diagonal, there is phase wrap- 
ping in some swings, e.g. clearly visible for i = 10°. 

One limitation of the rotating vector model is that it 
is difficult to distinguish between orientations with simi- 
lar |d(PA)/d/|. This occurs, for example, at two orientations 
with equal a and i — a. For a pure dipole at r <C ^lc, the 
Stokes phase portraits face a similar problem for a > 40°. 
Comparing (a, i) = (50°, 40°) and (50°, 60°), the hockey 
sticks in the I-Q plane are look alike. The balloons and 
heart shapes in the I-U and Q-U planes are slightly differ- 
ent, although the differences are close to / = and may 
be hidden if the real data is noisy. However, for a < 30°, 
the Stokes phase portraits clearly discriminate between dif- 
ferent orientations with equal i — a, in a way that the PA 
swing does not. For example, let us compare the orienta- 
tions (a,i) = (30°, 20°) and (30°, 40°) in Figure fTUl In the 
I-Q plane, the hockey stick at (a, i) = (30°, 20°) is more 
curved than at (a, i) = (30°, 40°). In the I-U plane, there is 
a distinct figure-eight instead of a balloon, whereas in the Q- 
U plane, there are two interlocking ovals instead of a heart 
shape. In this case, the Stokes phase portraits provide an 
obvious way to differentiate between the two orientations. 
We note that in Sections [4] and [5] the differences in the 
phase portraits for orientations with equal |i — q] are even 
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-0.5 0-0 0.5 1.0 1.5 --.0 -0.5 0.0 0.5 1.0 

-Q 




-0.5 0-0 0.5 1.0 1.5 -0.8 -0.2 0.2 0.8 1.2 



Figure 10. Two orientations with equal PA-swing gradients 
(|d(PA)/(ii|) at the inflection point, but different Stokes phase 
portraits. Top four panels: (a, i) = (30°, 20°), bottom four pan- 
els: (a,i) = (30°, 40°). Within each half of the figure, the four 
subpanels display (clockwise from top left): (a) 1/Ims.x (bottom) 
and PA (top; in radians) as functions of pulse longitude, (b) I-Q, 
(c) Q-U, and (d) I-U. 

more marked. Hence the Stokes phase portraits are useful 
supplements to the PA swing in many instances. 

The Stokes phase portraits produced by L = / cos 6o 
are extremely similar to those produced by a model with 
L oc I. For the sake of brevity, we do not investigate the 
latter model here. 

3.1.2 L = Ismeo 

Examples of the / (solid curve) and L (dashed curve) profiles 
for {a, i) — (30°, 30°) are shown in FigurelTTI / is normalised 
by its peak value. The magnetic colatitude 6o{t) of the emis- 
sion point (dotted curve) is also shown. The following trends 
are observed, (i) Along the a = i diagonal, where 6o goes to 
zero at the centre of the pulse, the sin 6o dependence results 
in a double-peaked L profile. Away from the diagonal, L/I 
increases, reaching 0.8 at (a, i) = (90°, 10°). (ii) Away from 
the diagonal, L follows 7, as in Section 13.1.11 resulting in 
similar phase portraits. 

The Stokes phase portraits for a filled core with L/I = 
sin Oq are drawn in Figures I12H14I Comparing the Stokes 
phase portraits and PA swings to those of Section 13.1.11 we 
note the following differences, (i) In the I-Q phase portrait 
(Figure [12)) , along the a = i diagonal, instead of a straight 
line, we see a 'flipped' hockey stick which slopes in the op- 
posite direction to the hockey sticks for a ^ i. (ii) In the 
I-U phase portrait (Figure [13]), along the a = i diagonal. 




Pulse longitude, L 

Figure 11. Dipole field at r ^ '"LC- Example of a pulse profile 
for a filled core beam with linear polarization L = /sin So and 
(a,i) = (30°, 30°). Solid, dashed and dotted curves represent the 
total polarized intensity /, degree of linear polarization L, and 
emission point colatitude Oq. Pulse longitude I is measured in 
units of degrees. 

there are broad balloon shapes instead of narrow ovals, and 
they extend leftwards from / = 1 as opposed to rightwards 
from 7 = 0. (iii) In the Q-U plane (Figure [14]) , instead of 
ovals, the diagonal now contains heart shapes with an ex- 
tremely pronounced cusp, narrowing with increasing a — i. 
There is no noticeable difference between the PA swings in 
this model (Figure llSp and those of Section 13.1.11 This is 
another example of how the Stokes phase portraits reveal 
differences in the beam geometry and polarization pattern 
which are not apparent from the PA swings. 

In summary, a filled core beam produces single-peaked 
pulses, resulting in several unique shapes in the Stokes 
planes. In the I-Q plane, a hockey stick shape is typically 
seen. In the I-U plane, we see a balloon shape which twists 
into a figure-eight in the range a ^ 50°, i 5^ 20°. In the Q-U 
plane, we see heart shapes for all i < a, and balloons for 
i~a^ 10°. 



3.2 Hollow cone emission 

3.2.1 L^IcosOo 

Examples of the 7 (solid curve) and L (dashed curve) profiles 
for (a,i) = (70°, 20°) and (70°, 60°) are shown in Figure [T6l 
7 is normalised by its peak value. The magnetic colatitude 
6o{t) of the emission point (dotted curve) is also shown. The 
following trends are observed, (i) Emission from a hollow 
cone results in double-peaked pulse profiles in the range \i — 
a\ ^ 30°, where the path traced by xo(t) crosses the edge 
of the cone twice, (ii) For a < 45°, the central dip becomes 
shallower as i increases. The reverse is true for a > 45°. (iii) 
We find L/I ~ 1 for a = i, decreasing to « 0.6 at {a,i) = 
(10°, 90°). The pulses are generally « 50% broader than in 
Section [3. 11 even though a is the same in both instances. 

The Stokes phase portraits for a hollow cone with 
L/I = cos 6o are drawn in Figures [T7t419l The paths of xo(i) 
shown in Figure [5] still apply but they cross a different beam 
pattern with a bright band at 9o = p. Hence the Stokes phase 
portraits look different to the left and right of the a = i di- 
agonal. 

In the I-Q phase portraits (Figure I17|) . we note the 
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Figure 6. Dipolc field at a low emission altitude r <C ^lc- Look-up table of Stokes phase portraits in the I-Q plane for a filled core 
beam with degree of linear polarization L = / cos 60 1 where 80 is the emission point colatitude. The panels are organised in landscape 
mode, in order of increasing 10° ^ i ^ 90° (left— right) and 10° ^ a ^ 90° (top— bottom) in intervals of 10°. I is plotted on the horizontal 
axis and normalised by its peak value. Q is plotted on the vertical axis. 
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Figure 7. Dipole field at r ^lc- Layout as for Figure [6l but for 1-lJ {I on tiie iiorizontal axis). 
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Figure 8. Dipole field at r <g; tlc- Layout as for Figure [B] but for Q-U {Q on the iiorizontal axis). 
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Figure 9. Dipole field at r <C rLC- Layout as for Figure|6] but for position angle (on the vertical axis in landscape orientation, in units 
of radians) versus pulse longitude (on the horizontal axis, in units of 27r radians). 
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Figure 12. Dipolc field at r <C tlc- Look-up table of Stokes phase portraits in tiie I-Q plane for a filled core beam with degree of linear 
polarization L = /sin^oi where 9q is the emission point colatitude. The panels are organised in landscape mode, in order of increasing 
10° < j ^ 90° (left-right) and 10° a ^ 90° (top-bottom) in intervals of 10°. / is plotted on the horizontal axis and normalised by its 
peak value. Q is plotted on the vertical axis. 
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Figure 13. Dipole field at r <C r'LC- As for Figure [121 but for I-U (I on horizontal axis). 
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Figure 14. Dipole field at r <C »'lc- As for Figure [121 but for Q-i 
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Figure 15. Dipole field at r <g; rLC- Layout as for Figure [TTl but for position angle (on the vertical axis in landscape orientation, in 
units of radians) versus pulse longitude (on the horizontal axis, in units of 27r radians). 
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100 200 300 

Pulse longitude. I 

Figure 16. Dipole field at r <C »'lc- Examples of pulse profiles 
for a hollow cone with opening angle 25° and linear polarization 
L = I cos So for two different orientations (a, i) (top and bottom 
panels). Solid, dashed and dotted curves represent the total po- 
larized intensity I, degree of linear polarization L, and emission 
point colatitude Bo- Pulse longitude is measured in units of de- 
grees. Top: double-peaked pulse with {ct,i) = (70°, 60°); bottom: 
single-peaked pulse with {a,i) = (70°, 20°). 



following trends, (i) Where the pulse is double-peaked, the 
hockey sticks twist into a tilted 7 shape. Where the pulse 
is single-peaked, the patterns are similar to Figure (6] (ii) 
For a = i, the hockey stick twists into a 7-shape as a = i 
increases. 

In the I-U phase portraits (Figure [18}, we note the fol- 
lowing trends, (i) Where the pulse profile is double-peaked, 
the I-U portraits feature two interlocking loops forming a 
triangular shape with twisted corners. Where the pulse is 
single-peaked, the patterns are similar to those of Figure [T] 
(ii) For a — i, the pattern evolves from a balloon to a trian- 
gular shape with an increasingly deep central cusp as a = i 
increases. 

Additional twisting in the patterns is also seen in the 
Q-U plane (Figure [19)1 . although the patterns are still sym- 
metric about (7 = 0. (i) For a = i, the patterns look quali- 
tatively similar to those of Section 13.1.21 although they are 
^ 90% broader, covering —1 ^ U ^ 1. (ii) For a > i, we 
see heart shapes and nested ovals, whereas for a < i, we see 
simpler, convex balloons. 

The PA swings (Figure [20[) are identical to the previous 
cases, the only difference being that they are visible above a 
threshold intensity over a greater fraction of the pulse period 
due to the broader pulse profiles. 




Figure 21. As for Figure [T6l but with linear polarization L = 
I sin 6. 



3.2.2 L = Isin6o 

Examples of the / (solid curve) and L (dashed curve) profiles 
for {a,i) = (70°, 20°) and (70°, 60°) are shown in Figure [2T] 
We note the following trends, (i) L follows / in the same 
way as Figure [T^ However, L/I peaks at « 0.5 on the a = i 
diagonal, (ii) L/I increases as \a — i\ increases, reaching « 
0.8 for \a~i\= 80°. 

The I-Q, I-U and Q-U phase portraits are shown 
in Figures I22H24I The phase portraits are similar to Sec- 
tion [33TT| except for orientations that yield double-peaked 
pulses. Here, the patterns are less twisted. Compared to Sec- 
tion [33II] the lower peak value of L/I leads to a narrower 
range of U and Q. For example, in the Q-U plane, along 
the a = i diagonal, the heart shape for (a,i) = (10°, 10°) 
is approximately 80% smaller than in Figure [TS] along both 
the U and Q axes. The PA swings (Figure I25p are exactly 
identical to Figure [20] 

In summary, a hollow cone beam pattern produces 
Stokes phase portraits that differ significantly from those of 
the filled core beam, when the pulses are double-peaked (for 
a — i| ^ 30°). In the I-Q plane, double-peaked pulses trace 
out a tilted 7 shape, whereas single-peaked pulses trace out 
a hockey stick. In the I-U plane, double-peaked pulses pro- 
duce a triangle with twisted corners, whereas single-peaked 
pulses produce a balloon. In the Q-U plane, the double- 
peaked pulses produce a mix of heart shapes and nested 
ovals, whereas the single-peaked pulses produce heart shapes 
(for a > i) and balloons (for a <i diagonal). 
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Figure 17. Dipolc field at r <C tlc- Look-up tabic of Stokes phase portraits in the I-Q plane for a hollow cone with opening angle 25° 
with degree of linear polarization L = I cos 6a, where do is the emission point colatitude. The panels are organised in landscape mode, 
in order of increasing 10° ^5 « ^ 90° (left-right) and 10° ^ a C 90° (top— bottom) in intervals of 10°. / is plotted on the horizontal axis 
and normalised by its peak value. Q is plotted on the vertical axis. 
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Figure 18. Dipole field at r <g; tlc- Layout as for Figure [TTl but for I-U (/ on fiorizontal axis). 
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Figure 19. Dipole field at r <g; tlc- Layout as for Figure [TTl but for Q-U (Q on horizontal axis). 
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Figure 20. Dipole field at r <g; rLC- Layout as for Figure [TTl but for position angle (on the vertical axis in landscape orientation, in 
units of radians) versus pulse longitude (on the horizontal axis, in units of 27r radians). 
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Figure 22. Dipolc field at r <C tlc- Look-up tabic of Stokes phase portraits in the I-Q plane for a hollow cone with opening angle 25° 
with degree of linear polarization L = I sin Sq , where 6[) is the emission point colatitude. The panels are organised in landscape mode, 
in order of increasing 10° ^5 « ^ 90° (left-right) and 10° ^ a C 90° (top— bottom) in intervals of 10°. / is plotted on the horizontal axis 
and normalised by its peak value. Q is plotted on the vertical axis. 
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Figure 23. Dipolc field at r <g; tlc- Layout as for Figure [22l but for I-V {1 on horizontal axis). 
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Figure 24. Dipole field at r <g; tlc- Layout as for Figure [22] but for Q-U (Q on horizontal axis). 
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Figure 25. Dipole field at r <C tlc- Layout as for Figure [22] but for position angle (on the vertical axis in landscape orientation, 
units of radians) versus pulse longitude (on the horizontal axis, in units of 27v radians). 
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4 DIPOLE FIELD AT r = 0.1 tlc 

Relativistic aberration affects the observed Stokes parame- 
ters in two ways. Firstly, the change in the emission point 
due to relativistic beaming [Eq. shifts the pulse centroid 
and hence I it). Beaming also shifts the PA swing in the 
same direction, resulting in a relative phase shift of 4r/rLc 
between the PA inflection point and the centroid of the pulse 
profile. Note that as we are only interested in the relative 
phase, we ignore time-delay effects, which phase-shift both 
the pulse pr ofile and the PA swing by an additional — r/rLC 
(|Dvksll2008l '). 

In this section, we present look-up tables for a dipole 
field at an emission altitude r — O.lrLC, where aberration 
distorts the Stokes phase portraits significantly while leaving 
the PA swings largely unchanged. We investigate the same 
beam patterns as in Section [3l but for the sake of brevity, 
we present look-up tables for L/I = cos So only. For each 
beam pattern (core emission in Section [4.11 conal emission 
in Section r4.2p we construct three phase portraits {I-Q, I-U 
and Q-U) for each pair of angles (a, i) at a fixed emission 
altitude (r — 0.1 tlc) and tabulate the results in Figures 

4.1 Filled core beam 

In this section, we point out the major differences between 
the Stokes phase portraits for a filled core beam and L = 
/cos So at r = O.lrLC fFigures l26H28p and the corresponding 
non-aberrated situation at r <C tlc (Figures [SHI]) . 

In the I-Q plane, the hockey sticks and straight lines 
in the non-aberrated case (Figure [HI broaden into balloons 
and ovals when aberration becomes important (Figure [26)) . 
The broadening increases as the panels approach the a = i 
diagonal. 

In the I-U plane (Figure I27|) . the symmetry about 
[/ = is broken, and the balloons and ovals are tilted. For 
a > i (below the diagonal), the balloons tilt upwards (major 
axis has dU/dl > 0), whereas for a < i (above the diago- 
nal), the ovals tilt downwards (major axis has dU/dl < 0). 
Another difference is that for Iq — j| = 10°, the phase por- 
traits feature figure-eights instead of ovals. Along the a — i 
diagonal there is a distorted Z-shape, although we note that 
the sharp breaks in the pattern are numerical artifacts. 

In the Q-U plane (Figure [28]) , the heart shapes in the 
a > i region (below the diagonal) are rotated counter- 
clockwise. As the panels approach the a = i diagonal, the 
top ventricle of the heart shape grows larger than the bot- 
tom ventricle. In the a < i region (above the diagonal), the 
ovals are rotated clockwise and are tilted downwards (major 
axis has dU/dQ < 0). Along the a = i diagonal, there are 
distorted heart shapes instead of ovals. 

4.2 Hollow cone 

We now examine the major differences between the Stokes 
phase portraits for a hollow cone beam and L — I cos 6o 
at r = O.lrLC (Figures I29H31|) and the corresponding non- 
aberrated case (Figures I17l419p . 

Similar to the filled core beam, aberration broadens 
the shapes in the I-Q plane (Figures [29]). At the orienta- 
tions where the pulse profile is double-peaked, the 7-shapes 



at Q > 40° and i > 40° from Figure [17] are replaced by 
mosquito shapes. For a < 30° and i < 30° , the I-Q portraits 
maintain their general shape but are broader. An exception 
is {a,i) = (20°, 30°), which changes from a C-shape in the 
non-aberrated case to a figure-eight with aberration. 

In the I-U plane, at the orientations where the pulse 
profile is double-peaked, the trefoils, triangles, and mosquito 
shapes retain their general shape. However, they are dis- 
torted and asymmetric about ?7 = 0, as the individual com- 
ponents of the patterns (e.g. the 'wings' of the mosquitoes) 
tilt at different angles and have different widths. 

The shapes in the Q-U plane do change at those ori- 
entations where the pulse profile is double-peaked. In the 
a > i region (below the diagonal), the interlocking ovals 
become distorted, asymmetric heart shapes as a and i in- 
crease. In the a > i region (above the diagonal), there are 
figure-eights [e.g. {a,i) = (20°, 40°)] which become distorted 
heart shapes as a and i increase. Along the a = i diagonal, 
distorted heart shapes evolve into distorted mosquito shapes 
as a and i increase. 

4.3 How widespread are low emission altitudes? 

As our first application of the above theory to observa- 
tional data, we ask: what proportion of pulsars with super- 
ficially dipole-like characteristics (e.g. S-shaped PA swings) 
are truly consistent with a pure dipole field and low, non- 
aberrated emission altitude (Section [3]l , when we examine 
their Stokes phase portraits? 

To this end, we select two sets of supe rficially dipole- 
like pulsars from the 366 objects pre sented in lGould fc Lvnj 
(1998) and iManchester et al.l (| 19981 ). Set I co ntains pulsars 
that follow closely the W oc P~^/'^ relation (lRankinlll99"3h 
and hence are likely prima facie to have an approximately 
dipolar field in the emission region, with or without aberra- 
tion. Tlie^'_aje_djawn_^|om near the diagonal line in Figure 
9 of lGould fc Lvnd { 19981 ): as a rule, their pulse width at 
0.1/max (where /max is the peak intensity) is less than 20°. 
Set II contains pulsars that have a clean, S-shaped PA swing, 
after allowing for phase wrapping. From the suitable pulsars, 
we choose 12 at random in each set. 

We present the pulse profiles. Stokes phase portraits, 
and PA swings for Set I at 610 MHz in Figures [32] and 
1331 PA data points are only plotted for L ^ O.lLmax and 
/ ^ O.l /max- The data are sourced from the EPN online 
archive l|Lorimer et al.lll998l ). The 12 objects are calibrated 
separately, so the polarization basis varies from pulsar to 
pulsar. For the purposes of this section, we do not attempt 
to determine /3 for each object; we focus on the shape of the 
Q-U patterns, not their orientations. From Sections [3] and 
13] we note that one clear fingerprint of low-altitude emission 
from a pure dipole is the symmetry of the Q-U pattern about 
U = Q (rotated in general by an angle 2/3) . Hence, in Figures 
I32H35I any pattern which is not symmetric about any line 
either (i) cannot originate from a low altitude where aberra- 
tion is negligible, or (ii) is the result of a largely tilted axis 
of symmetry, as discussed in Section r2.4.1l It is not possible 
to disentangle these effects solely from linear polarization 
data. Circular polarization phase portraits may assist, but 
their analysis lies outside the scope of this paper. 

Of the 12 objects in Figures l32l and [33l four have Stokes 
phase portraits which roughly resemble a dipole at low alti- 
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Figure 26. Dipolc field at r = O.lrLC- Look-up tabic of Stokes phase portraits in the I-Q plane for a filled core beam with degree 
of linear polarization L = /cos 60, where Oq is the emission point colatitude. The panels are organised in landscape mode, in order 
of increasing 10° ^ « ^ 90° (left— right) and 10° C a ^ 90° (top— bottom) in intervals of 10°. / is plotted on the horizontal axis and 
normalised by its peak value. Q is plotted on the vertical axis. 



28 C. T. Y. Chung et al. 



A 



6 




] >V1 U □> I 



A 



A 



A 



A 



A. 



A. 



A 



A. 



A 



A 



A 



A. 



A 



A 



/I 



^ 



^ ^ S ' 



O p p ■ 



U D □ O O 



■5 ■& ^ & 



pop 
e ki ^» 



A 



' e- p p 



" h P P P 

I- ^ 1.; I4J ^ 



,pfhpp 



.-1 



V 



■5 !& * & * 



n 4^. & I ^ tn W ^ ^ 



./J 



U'' 



■ O O O D 



/ / 



o o o □ o 



A 



A 



Figure 27. Dipole field at r = O-Itlc- Layout as for Figure [26l but for I-U (1 on the horizontal axis). 
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Figure 28. Dipolc field at r = O.Itlc- Layout as for Figure [26l but for Q-U {Q on tiie horizontal axis). 
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Figure 29. Dipolc field at r = O.lrLC- Look-up tabic of Stokes phase portraits in the I-Q plane for a hollow cone with opening angle 25° 
with degree of linear polarization L = / cos 6*0, where 9a is the emission point colatitude. The panels are organised in landscape mode, 
in order of increasing 10° ^ i ^ 90° (left-right) and 10° ^ a C 90° (top— bottom) in intervals of 10°. / is plotted on the horizontal axis 
and normalised by its peak value. Q is plotted on the vertical axis. 




Figure 30. Dipole field at r = O.lrLC- Layout as for Figure [29l but for I-U (/ on the horizontal axis). 
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Figure 31. Dipole field at r = O.lrLC- Layout as for Figure [29l but for Q-U (Q on tiie iiorizontal axis). 
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tudes. PSR J0629+2415 and PSR J1740-3015 (Figure [32l 
third and fifth rows) match a core beam for large i and small 
Q. Both L = I cosOo and L = I sin 80 give similar phase por- 
traits for these angles. PSR J1823-3106 and J1926+1648 
(Figure [33l first and fourth rows) match a core beam with 
a ^ i and L — Jcos^o- The other eight pulsars have Q- 
U portraits which are significantly asymmetric. This is re- 
flected in their PA swings, which are not as clean as the other 
four pulsars. For example, PSR J2046-K1540 (Figure [551 fifth 
row) features an asymmetric mosquito in the Q-U plane. As- 
suming that the asymmetry is the result of relativistic aber- 
ration, such a pattern cannot be reconciled with emission 
from low altitudes. It matches a dipole field at r = O.lrLC 
along a = i for a hollow cone with L = /cos So (see the 
look-up table in Figure [5T|) . where aberration is important. 
Alternatively, if the asymmetry is the result of a tilted axis of 
beam symmetry, it is possible that the emission originates 
instead from a low altitude, and that aberration does not 
play a part in distorting the phase portraits. 

Figures [34] and [35] display the pulse profiles, Stokes 
phase portraits, and PA swings for Set II. This set com- 
prises data at raultiple frequencies (|Gould fc Lvnel 1 19981 : 
[Manchester et aU Il998h . sourced from the EPN online 
archive. Four of the 12 pulsars in this set feature approx- 
imately symmetric mosquito and heart shapes matching a 
dipole at low emission altitudes. PSR J0528-f2200 (Figure 
\M second row), PSR J2048-1616, and J2346-0609 (Figure 
1351 third and sixth rows) match the Q-U portraits for a hol- 
low cone with a = i ^ 30° {L = I cos 9o and L = I sin 60 give 
similar answers for these angles). PSR J2113+4644 (Figure 
1351 fourth row) matches all three phase portraits for a filled 
core beam for 0° ^ a ~ i ^ 10° (again, L = I cos 9q and 
L = I sin 60 give similar answers for these angles) . In Figure 
1341 and [35l the mosquito and heart shapes of the other eight 
pulsars are noticeably asymmetric. Again, we note that the 
asymmetry may be the result of a tilted axis of beam sym- 
metry, or it may indicate that the emission originates from 
higher altitudes, where aberration takes effect. 

For a dipole field (or indeed any axisymmetric field), 
aberration changes the relative phase between the PA swing 
and pulse profile, but it does not cause the shape of the 
PA swing to change with emission frequency (unlike the 
pulse width). In reality, multi- frequency ob servations show 
that P A swings do change with frequency l|johnston et al.l 
l2008bl ). We return to this issue in Section [5.3l where we find 
that the Stokes phase portraits for magnetic geometries with 
an appreciable toroidal component at r > O.lrLC evolve with 
frequency dramatically and in an informative way. 

A key tenet of pulsar radio emission theory holds that 
core and cone emission mostly comes from low altitudes, 
where departures from a dipole field are small. The results 
in Figures [321351 implv that in a clear majority (> 60%) of 
objects, the emission may well originate from high altitudes, 
where aberration and/or the toroidal field component are 
significant. We now turn to explore this important issue in 
Section [5] 

5 CURRENT-MODIFIED DIPOLE FIELD 

In this section, we examine the pulse profiles. Stokes phase 
portraits, and PA swings associated with a current-modified 



magnetic field composed of a dipole with symmetry axis 63 
plus a toroidal component (cylindrically symmetric about 
63) with magnitude 

-B^ = —-Bp cos Q sin Sr/rLC, (16) 

where Bp = (B^ + Bg)^^^ is the poloidal field strength. 
Equation (fT6)l is the scaling expected for r < tlc if B^ 
is generated by a field-aligne d Goldreich-Julian current in a 
plasma- filled magnetosphere (|Hibschman fc Arons|[200ll ). If 
B^ is generated by the displacement current, as in a vacuum 
rotator, it scales as r^; we save this case for future analysis 
(jPe utsc h 1955; Melato§i997.). 

We investigate the same beam and polarization patterns 
as in Section [3] The emission point xo(f) traces out a differ- 
ent locus in the current-modified dipole as compared to the 
pure dipole, and the locus varies differently with altitude. In 
Sections l5.1l and l5.2l we present two-dimensional look-up ta- 
bles of Stokes phase portraits for orientations 10° 5S i ^ 90° 
and 10° ^ Q ^ 90°. For each beam pattern (core emission 
in Section 15.11 conal emission in Section 15. 2p and polariza- 
tion pattern {L/I — cos So, sin So), we construct three phase 
portraits (I-Q, I-U and Q-U) for each pair of angles {a, i) 
at a fixed emission altitude (r — 0.1 tlc). In Section [5.31 we 
show some examples of how the phase portraits change with 
altitude. 



5.1 Filled core emission 

5.1.1 L = /cosSo 

An example of the I (solid curve) and L (dashed curve) 
profiles for a single-peaked pulse with {a, i) = (70°, 30°) are 
shown in Figure [551 / is normalised by its peak value. The 
dotted curve shows how the magnetic colatitude 9o{t) of the 
emission point varies across one pulse period. We note the 
following trends, (i) The emission points and pulse profiles 
are phase-shifted by approximately — r/rLC radians due to 
aberration, (ii) The profiles behave similarly to the pure 
dipole (Section I3.1.ip . narrowing with increasing a and i. 
(iii) For |q! — i| ^ 10°, 0o remains near zero for some time 
0.5 phase units at (q, i) — (10°, 10°) and ~ 0.13 phase 
units at {a,i) = (90°, 90°)]. This is a numerical artifact 
caused by the finite grid resolution near So = 0. 

Figure [37] shows the path traced by xo(i) in the body 
frame across one pulse period for three cases: (A) (a, i) — 
(20°, 10°), (B) (20°, 30°), and (C) (20°, 50°). As in Figure 
[5] the path changes from an undulation for a < i (curve 
C) to an ellipse for a > i (curve A). Note that curve A 
features a cusp near {0,4>) « (0.4 rad, 1.6 rad). The toroidal 
magnetic field twists the paths, so that they are not refiec- 
tion symmetric about some longitude, unlike in Figure [S] 
The asymmetry distorts the Stokes phase portraits, but the 
clear distinction between a < i and a > i persists. 

The Stokes phase portraits are drawn in Figures [5511401 
The phase portraits differ from those of the pure dipole in 
both Sections [3] and [4] In the I-Q plane (Figure [38}, we note 
the following behaviour, (i) For a — i, the pulse traces out 
a triangular shape. The cutoff of the shapes at 7 = 1 is a 
numerical artifact caused by the finite grid resolution near 
So = 0. (ii) For a > i, we see tilted, narrow ovals. The ovals 
tilt upwards (major axis has slope dQ/dl > 0). For i ^ 20° 
and a > 45°, the balloon twists into a figure-eight, (iii) For 
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Figure 32. Pulse profiles, Stokes phase portraits and PA sw ings for a selection of six superficially dipolar pulsars that closely follow the 
pulse-width-period relation at 610 MHz llGould &: Lvnelll998^ . The data for each pulsar occupy a row each in landscape mode. From left 
to right, the columns show: (1) / (mjy, solid curve) and L (mjy, dashed curve), (2) the I-Q phase portrait, (3) the I-U phase portrait, 
(4) the Q-U phase portrait, and (5) the PA swing (data points with L > O.lLmax and / ^ 0.1/max plotted only). Data arc presented 
courtesy of the EPN online archive. 
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Figure 33. As for Figure [32l for six other objects selected according to the same criterion. 
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Figure 34. Pulse profiles, Stokes phase portraits and PA sw ings for a selection of six superficially dipolar pulsars that exhibit a 
clean, S-shaped PA swi ng. All pulsars are obser ved at 610 MHz llGould fc Lvnelll998h . except for PSR J0536-7543 (663 MHz) and PSR 
J1751-4657 (433 MHz) ^Manchester et al.|[T99l) . The data for each pulsar occupy a row each in landscape mode. From left to right, the 
columns show: (1) / (mjy, solid curve) and L (mjy, dashed curve), (2) the I-Q phase portrait, (3) the I-U phase portrait, (4) the Q-U 
phase portrait, and (5) the PA swing (data points with L ^ O.lLmax and / ^ O.l/max plotted only). Data are presented courtesy of the 
EPN online archive. 
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Figure 35. As for Figure 1321 for six other objects selected according to the sam e criterion. All pulsars are observed at 610 MHz 
llGould fc Lvnelll998l'). except for PS R J1842-0359 (1.408 GHz) l lGould fc LvnJlQQ^'l . PSR J2144-3933 (659 MHz) and PSR J2346-0609 
(435 MHz') l lManchester et al.lll998h . 
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a < i, we see tilted balloons for |a — i| ^ 45°, which narrow 
into tilted ovals as i increases. The balloons tilt downwards 
(major axis has dQ/dl < 0) for i ^ oi < 25° and upwards 
(major axis has dQ/dl > 0) for i — a > 25°. 

In the I-U plane (Figure [39)) . we note the following be- 
haviour, (i) For a = j, we obtain a triangular shape, (ii) For 
a > i (below the diagonal in Figure [39)) . the pulse traces out 
a balloon which narrows into a straight line as a increases, 
(iii) For a < i (above the diagonal), we obtain balloons for 
i — Q < 10° , and narrow ovals for i — a > 20° . 

In the Q-U plane (Figure I40p . we see a mix of bal- 
loons, heart shapes and ovals with the following proper- 
ties, (i) Along the a = i diagonal, we observe a 'C shape 
for (a, i) = (10°, 10°), which evolves into a distorted heart 
shape as a = i increases, (ii) For a > i (below the diago- 
nal in Figure [40)) and |ia — i| < 25°, we see distorted heart 
shapes with a cusp at (Q,U) = (0,0). For |q — i| > 25°, 
we see tilted balloons which narrow into a straight line as a 
increases, (iii) For a < i (above the diagonal), we see heart 
shapes for |a — j| ^ 10°. These shapes evolve into tilted bal- 
loons as |a — i| increases, (iv) As in Figure [5^ the patterns 
are broadest for |a — i| ^ 10° and narrow as Iq — i| increases. 

Generally, the shapes in Figure |40] are not symmetric 
about U = or Q = and have a higher degree of rotation 
than those in Section |4| The amount by which the shapes 
rotate depend on |q — i|. Geometrically, this happens because 
6o{t) increases with |a — i|. The degree of rotation decreases 
as |a — «| increases. 

In Figure 1411 we plot the PA swings corresponding to 
each panel in Figures I38H40I We plot only the parts of the 
swing that are illuminated by the pulse, i.e. when L 10~^. 
We note the following behaviour, (i) There are clear distor- 
tions in the S-shape for \a — i\ ^ 10°. Along a = i, where 
the I ~ 1 saturation occurs, d{PA)/dl changes discontin- 
uously. A good example is {a,i) = (20°, 20°). The kinks 
around the discontinuities are numerical artifacts, but the 
smoothly varying curves leading up to the kinks are real, 
(ii) For |a — i| ^ 20°, the shape of the PA swings appear 
similar to those of Section^ The toroidal field enhances the 
phase shift of the PA swings in Figure [41] For example, at 
(a, i) = (40°, 70°), the relative phase shift between the PA 
swing and the pulse profile is « 5r/rLC- The phase shift is 
altitude dependent; we discuss altitudinal variations in Sec- 
tion 15.31 The similarity of the PA swings for the pure and 
current-modified dipoles reinforces the necessity to supple- 
ment PA swings with Stokes phase portraits when diagnos- 
ing the pulsar magnetosphere. 

5.1.2 L = Ismeo 

An example of the I (solid curve) and L (dashed curve) pro- 
files for (q, i) = (30°, 30°) is shown in Figure[42l The profiles 
are similar to those of the pure dipole in Section 13.1.21 dis- 
playing the same double-peaked L profile along the a = i 
diagonal. The only difference is that the path traced by the 
emission point is not symmetric about the pulse centroid. 

The Stokes phase portraits for a filled core with L/I = 
sin^o are drawn in Figures 14314451 Comparing the Stokes 
phase portraits and PA swings to those in Section 15.1.11 we 
note the following differences, (i) In the I-Q phase portrait 
(Figure I43|) . along the a — i diagonal, there is a balloon 
instead of a triangle. The panels where L/I is smaller than 




Pulse longitude, L 

Figure 36. Current-modified dipole field at r = O.lrLc- Example 
of a pulse profile for a filled core beam with linear polarization 
L = /cosfio and {a,i) = (70°, 30°). Solid, dashed and dotted 
curves represent the total polarized intensity I, degree of linear 
polarization L, and emission point colatitudc Oq. Pulse longitude 
I is measured in units of degrees. 
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Figure 37. Current-modified dipole field at r = O.lr'LC- Ex- 
amples of the path traced by the emission point xo(t) in the 
body frame (curves), overplotted on the intensity map for a filled 
core beam (greyscale). The beam is brightest at 9 = (cor- 
responding to the 63 axis). Curves A: (o,i) = (20°, 10°), B: 
ia,i) = (20°, 30°), and C: (o,i) = (20°, 50°). 

in Section 15.1.11 exhibit narrower shapes than their coun- 
terparts in Figure 1381 for example, along the diagonal, the 
shapes are ~ 10 times narrower than in Figure 1381 (ii) In 
the I-U phase portrait, (Figure [44)) . the 1 = 1 wall seen in 
Figure [35] for |q — i| ^ 10° is no longer present. For a = i, 
instead of triangles, we see balloons. In the adjacent pan- 
els, we see tilted balloons which twist into a figure-eight for 
a 20° and a ^ 60°. Again, the shapes along the a = i 
diagonal are « 10 times narrower than their counterparts in 
Figure [39] (iii) In the Q-U plane (Figure [45j , the main dif- 
ference between Figures [ID] and [^ occurs for \a — i\ = 10°, 
where the heart shapes are distorted. Along the a = i diag- 
onal, we see a distorted figure-eight. There is no noticeable 
difference between the PA swings in this model (Figure [JSJ 
and the model in Section [5. 1.1 1 

In summary, the Stokes phase portraits for a filled 
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Figure 38. Current-modified dipole field at r = O.lrLC- Loolc-up table of Stolces pliase portraits in tlie I-Q plane for a filled core beam 
with degree of linear polarization L = / cos So, where do is the emission point colatitude. The panels are organised in landscape mode, 
in order of increasing 10° C * C 90° (left-right) and 10° $5 a ^ 90° (top-bottom) in intervals of 10°. I is plotted on the horizontal axis 
and normalised by its peak value. Q is plotted on the vertical axis. 
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Figure 39. As for Figure l38l but for I-U {I on the horizontal axis) 



V. 



J 



k 

i.. V 



J 



/ 



/ 



/ 



i.. v 



/ 



J 



/ 



.4 



J 



J 



/ 



/ 



mi 



/ 



Stokes tomography of radio pulsar magneto spheres. I. Linear polarization 



K 



'■—.J 



^ 



^ ^ s ' 



'rh p p., 

I 15 S "J" < 



2 o o c- a 



in 



n 



p ^ ^ M g 



n 



15 a Jii 



.<7 



^ gi J. < 



I □ o o - 



..■■-■~r 



// 



Figure 40. As for Figure [38l but for Q-U (Q on the horizontal axis). 
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Figure 41. Current-modified dipole field at r = O.lrLC- Layout as for Figure [38l but for position angle (on the vertical axis in landscape 
orientation, in units of radians) versus pulse longitude (on the horizontal axis, in units of 27r radians). 
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Pulse longitude, i 

Figure 42. Current-modified dipole field at r = O.lrLC- Example 
of a pulse profile for a filled core beam with linear polarization 
L = /sin So and (a,i) = (30°, 30°). Solid, dashed and dotted 
curves represent the total polarized intensity /, degree of linear 
polarization L, and emission point colatitude da- Pulse longitude 
I is measured in units of degrees. 



core beam are significantly different for pure and current- 
modified dipoles. Depending on its orientation, the toroidal 
field can either enhance or diminish the effects of aberra- 
tion. The oval, balloon, and heart shapes are not symmetric 
about any axis and are rotated about {U,Q) — (0,0). The 
PA swings are distorted for |q — i| ^ 10° but are similar to 
the pure dipole for other orientations apart from the phase 
shift. 



5.2 Hollow cone emission 

5.2.1 L = Jcos6lo 

Examples of the / (solid curve) and L (dashed curve) profiles 
for (q, i) = (70°, 30°) and (70°, 60°) are shown in Figure[47l 
I is normalised by its peak value. The magnetic colatitude 
do{t) of the emission point (dotted curve) is also shown. 
The profiles are double-peaked for |a — i| ^ 30° and behave 
similarly to the pure dipole at low altitudes fSection l3.2.ip . 
with the exception of the —r/r-LC phase shift. 

The Stokes phase portraits for a hollow cone with 
L/I = cos So are drawn in Figures [48U50I In the I-Q phase 
portraits (Figure [48]), we note the following trends, (i) For 
a = z, we obtain a tilted crescent at (a,j) = (10°, 10°), 
which evolves into an asymmetric mosquito shape as a — i 
increases, (ii) For a > i (below the diagonal in Figure [48)) . 
where the pulse is double-peaked, we obtain a mix of figure- 
eights and distorted mosquito shapes, (iii) For a < i (above 
the diagonal), where the pulse is double-peaked, the pulse 
traces mosquito shapes (for a > 25°) and a mix of balloons, 
figure-eights and triangular shapes (for a < 25°). 

In the I-U phase portraits (Figure [49]), we note the 
following trends, (i) Along the a = i diagonal, we obtain 
a tilted crescent at {a,i) = (10°, 10°) which twists into a 
mosquito shape as a — i increases, (ii) For a > i (below 
the diagonal), where the pulse is double-peaked, we obtain 
heavily distorted, complex shapes at a ^ 30°, which evolve 
into a distorted mosquito shape as a increases, (iii) For a < i 
(above the diagonal), where the pulse is double-peaked, we 




Pulse longitude, I 

Figure 47. Current- modified dipole field at r = O.lrLc. Exam- 
ples of pulse profiles for a hollow cone with opening angle 25° 
and degree of linear polarization L = IcosOq. Solid, dashed and 
dotted curves represent the total polarized intensity /, degree 
of linear polarization L, and emission point colatitude Oq. Pulse 
longitude I is measured in units of degrees. Top: double-peaked 
pulse with (o, i) = (70°, 60°); bottom: single-peaked pulse with 
(a, j) = (70°, 30°). 

obtain a tilted figure-eight at a = 10, which evolves into 
twisted, interlocking ovals as a increases. 

In the Q-U phase portraits (Figure I50|) . we note the 
following trends, (i) Along the a = i diagonal, the pattern 
evolves from a tilted heart shape to a figure-eight as q = i 
increases, (ii) Interestingly, in the panels to the left of the 
a = i diagonal, one of the ventricles of the heart shape 
twists into a secondary loop which rotates counter-clockwise 
as a increases [compare [ot,i) = (40°, 30°), (60°, 50°), and 
(90°, 80°)]. (iii) For a < i, where the pulses are double- 
peaked, a tilted figure-eight pattern emerges, (iv) The phase 
portraits for orientations which produce single-peaked pulses 
behave like in Section [5TT71] i.e. they rotate about {U,Q) = 
(0,0) as \a — i\ increases. 

The PA swings (Figure [ST]) are identical to the previous 
cases for all orientations, the only difference being that they 
are visible above a threshold intensity over a greater fraction 
of the pulse period due to the broader pulse profiles. 

5.2.2 L^IsinOa 

Examples of the I (solid curve) and L (dashed curve) profiles 
for (q, i) = (70°, 30°) and (70°, 60°) are shown in Figure [52] 
The profiles are similar to those of Section [3.2.21 
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Figure 43. Current-modified dipole field at r = O.lrLC- Loolc-up table of Stolces pliase portraits in tlie I-Q plane for a filled core beam 
and degree of linear polarization L = /sin^o, where Bq is the emission point colatitude. The panels are organised in landscape mode, in 
order of increasing 10° ^ i ^ 90° (left-right) and 10° ^ Q C 90° (top-bottom) in intervals of 10°. / is plotted on the horizontal axis and 
normalised by its peak value. Q is plotted on the vertical axis. 
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Figure 44. As for Figure B3l but for I-U {I on the horizontal axis). 
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Figure 45. As for Figure l43l but for Q-U {Q on the horizontal axis). 
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Figure 46. Current-modified dipole field at r = O.lrLC- Layout as for Figure l43l but for position angle (on the vertical axis in landscape 
orientation, in units of radians) versus pulse longitude (on the horizontal axis, in units of 2-k radians). 
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Figure 48. Current-modified dipole field at r = 0.1 r'LC- Look-up table of Stokes phase portraits in the I-Q plane for a hollow cone with 
opening angle 25° and degree of linear polarization L = /cos^Oi where So is the emission point colatitude. The panels are organised in 
landscape mode, in order of increasing 10° i 90° (left-right) and 10° ^5 a ^ 90° (top-bottom) in intervals of 10°. I is plotted on the 
horizontal axis and normalised by its peak value. Q is plotted on the vertical axis. 
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Figure 49. As for Figure l48l but for I-U {I on tfie horizontal axis). 
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Figure 50. As for Figure l48l but for Q-U (Q on the horizontal axis). 



Stokes tomography of radio pulsar magneto spheres. I. Linear polarization 51 



/ 



I 



1 



■7 



■7 



/ 



7 



T 



7 



o Qp o Q Q go a PC o o □ pp 



r 



O PP P DP PP 



o PP PP P PP o PPP P P PP o PP P P P DC 



7 



/ 



7 



ti § bttt^ Q bbbba a^bblsbb^ a^bbSbb^ Q^btziobbb QCiblbbbb^ 



/ 



7 



bb b Q Q PP p -abb abba 







■fi! is& jab 



o^t?^3bbbb D b^ b b b b b q i 



"A 



^ bo b Q b bp p bb b p Q bo p abb abba p ob q o b b Q 



obbbbbbb obbbbbbb c-bbc 



Figure 51. Current-modified dipole field at r = 0.1 ^lc- Layout as for Figure [48] but for position angle (on the vertical axis in landscape 
orientation, in units of radians) versus pulse longitude (on the horizontal axis, in units of 27r radians). 
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Pulse longitude, L 

Figure 52. Current-modified dipole field at r = O-Itlc- Exam- 
ples of pulse profiles for a hollow cone with opening angle 25° 
and degree of linear polarization L = I svciOq. Solid, dashed and 
dotted curves represent the total polarized intensity I, degree 
of linear polarization L, and emission point colatitude Oq. Pulse 
longitude I is measured in units of degrees. Top: double-peaked 
pulse with {a,i) = (70°, 60°); bottom: single-peaked pulse with 
(«, j) = (70°, 30°). 

The I-Q, I-U, and Q-U phase portraits are shown 
in Figures I53I455I All three phase portraits are similar to 
Section 15.2.11 except for orientations that produce double- 
peaked pulses. Here, the patterns are less twisted. The lower 
peak value of L// as compared to Section [5.2. II is accompa- 
nied by a narrower range of U and Q. For example, in the 
Q-U plane (Figure [55]) for (a, i) = (10°, 10°), the oval has 
approximately 20% the size of its counterpart in Figure [50l 
along both the U and Q axes. 

The PA swings in Figure are identical to Figure [5T1 
except in the squared-ofF valleys where 6o is near zero and 
we do not plot the curve {L < 10~^). 

5.3 Emission altitude 

In this section, we investigate how the emission altitude 
changes the Stokes phase portraits and PA swings. For the 
sake of definiteness, we focus on the orientation (a, i) = 
(70°, 30°) for a filled core beam with polarization L = 
7 cos 00- In Figure [57l we present the Stokes phase portraits, 
PA swing, and locus of xo(i) for this orientation at four al- 
titudes ranging from rmin ~ 0.05 tlc to rmax ~ 0.3 tlc- The 
pure dipole case at r ^ tlc is presented at the top of the 
figure as a comparison. / is normalised by its peak value. 



We begin with an obvious yet fundamental observation: 
in a current-modified magnetic dipole field, the phase por- 
traits evolve with altitude. The PA swing evolves also but 
to a lesser degree. 

We note the following trends, (i) The peak value of L/I 
increases from 0.90 at rmin to 0.92 at rmax. (ii) In the I- 
Q plane, at r = rmin, the hockey stick straightens into a 
line (major axis has dQ/dl > 0), ranging over < Q < 
0.8. As r increases, the slope of the line becomes less steep. 
At r = rmax, the line broadens into a tilted banana shape 
ranging over —0.07 < Q < 0.15. (iii) In the I-U plane, the 
balloon narrows and tilts upwards (major axis has dU /dl > 
0). (iv) In the Q-U plane, the heart shape narrows into a 
straight line, then broadens again as r increases, rotating 
counter-clockwise about {Q,U) — (0,0). (v) The shape of 
the PA swing changes as r increases, and the pulse longitude 
where the phase wrapping occurs shifts gradually from ~ 
0.65 phase units at r = rmin to ~ 0.54 phase units at r = 

rmax . 

To explain the trends, we look at the locus of xo(i) in 
the body frame. From r = rmin to r = rmax, xo(t) traces out 
an oval containing a secondary loop, which is not present 
for r << tlc- The oval tilts and becomes egg-shaped with 
increasing r. The range of 6o{t) increases from 0.4 rad < So 
1.3 rad at r = rmin to 0.4 rad < So ^ 1.4 rad at r = rmax. 
00 (i) shifts from 0.8 rad < (/)o ^ 2.2 rad at r = rmin to 
0.6 rad < 0o < 2.1 rad at r = rmax- As I and L do not 
depend on 00, the change in the shapes must be caused by 
the tilt and asymmetry of the locus. For a large tilt, the 
path traced from Z = 0to/ = 0.5 traverses different values 
of 9o than the path traced from Z = 0.5 to Z = 1, resulting in 
different ranges of and |?7| in the first and second halves 
of the pulse. This causes the asymmetry about the (/-axis 
and the broadening and narrowing of the shapes in the I-Q, 
I-U, and Q-U planes. 

Figure [57] refers to a particular orientation and polar- 
ization model, but the trends it depicts are fairly generic. 
We find that overall, from rmin to rmax, the phase portraits 
change modestly with emission altitude at a fixed {a,i). 
Their sizes, orientation, and detailed substructure (e.g. extra 
twisting of the secondary loops) change gradually, but there 
are no cases of the phase portraits changing drastically, for 
example from a hockey stick to a trefoil. Therefore, when 
comparing the look-up tables to real data, we do not have 
to worry about very degenerate matches, i.e. equally good 
fits for very diflerent combinations of (a, i) at two or more 
very different altitudes. We do, however, have to be aware 
that similar shapes like the hockey stick and banana can 
interchange as the altitude varies, leading to mild (r, a, i) 
degeneracies. 



6 CONCLUSION 

In this paper, we introduce Stokes phase portraits as an 
observational probe of a radio pulsar's magnetic geometry, 
orientation, and emission altitude. Our aim is to supply 
the reader with a set of easy-to-use tools for analysing pul- 
sar polarization data in a physically instructive way. We 
present handy look-up tables of pulse and polarization pro- 
files. Stokes phase portraits, and PA swings as functions 
of magnetic inclination angle a and line-of-sight inclination 
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Figure 53. Current-modified dipole field at r = 0.1 r'LC- Look-up tabic of Stokes phase portraits in the I-Q plane for a hollow eone with 
opening angle 25° and degree of linear polarization L = /sin^Oi where 9q is the emission point colatitude. The panels are organised in 
landscape mode, in order of increasing 10° S5 i ^ 90° (left-right) and 10° a ^ 90° (top-bottom) in intervals of 10°. I is plotted on the 
horizontal axis and normalised by its peak value. Q is plotted on the vertical axis. 
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Figure 54. As for Figure [53l but for I-U {I on the horizontal axis). 




Figure 55. As for Figure [53] but for Q-U {Q on the horizontal axis). 
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Figure 56. Current-modified dipole field at r = 0.1 ^lc- Layout as for Figure [53] but for position angle (on the vertical axis in landscape 
orientation, in units of radians) versus pulse longitude (on the horizontal axis, in units of 2-k radians). 
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Figure 57. Stokes phase portraits, PA swings and emission point locus xo(t) in the body frame for a filled core beam with polarization 
L = / cos So and = (70°, 30°), at five altitudes. In landscape mode, the top row shows a pure dipole at r <S ''lc- rows 2—5, r 

ranges from r = 0.05 tlc (second row) to r = 0.3 tlc (bottom row). From left to right in landscape mode, the columns show: (1) the 
I-Q phase portrait, (2) the I-U phase portrait, (3) the Q-U phase portrait, (4) the PA swing (data points with L ^ 10"'^ plotted only), 
and (5) the locus of xo(t). The intensity map (greyscale) follows Figure [51 
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angle i, for a pure and current-modified dipole, emission 
from a filled core and hollow cone, and two simple models 
of the degree of linear polarization. For the current-modified 
dipole, we also investigate how the phase portraits change 
with emission altitude r. In all the models, relativistic aber- 
ration is incorporated fully to order C(r/r-Lc)- To assist the 
reader, we summarise in Table [T] where to find the look-up 
tables corresponding to each model. 

For a pure dipole at low emission altitudes (r <C ^lc), 
the Stokes phase portraits display relatively simple shapes, 
summarised in Table [51 The shapes are classified according 
to colloquial descriptors in the Appendix. They are sym- 
metric about U = Q, because the effects of aberration are 
negligible. The curvature in the shapes is determined by the 
path traced by xo(i); for example, as \a — i\ increases, the 
curvature of the hockey stick in the I-Q plane decreases. For 
a filled core, the phase portraits produced by the two polar- 
ization models, L = I cos do and L = / sin So are generally 
similar, except for oi = i, where L is double-peaked. When 
the filled core is replaced by a hollow cone, the pulse pro- 
files for |q — i| ^ 30° become double-peaked, and the phase 
portraits transform into more complicated shapes. 

For a pure dipole at r = O.lrLC, aberration distorts the 
Stokes phase portraits. The symmetry about ?7 = is bro- 
ken, causing the shapes in the I-U and Q-U planes to tilt and 
rotate. In the I-Q plane, the shapes broaden noticeably. A 
phase shift between the pulse centroid an d PA swing infiec- 
tion point of » 4r/rLC is also observed (|Blaskiewicz et al.l 
Il99ll :[ Hibschman fc Aronj|200ll : lDvksll2008ll . 

For a current-modified dipole at r = O.lrLC, the phase 
portraits are also asymmetric about Q = Q and U = Q. 
The toroidal field can either partially cancel out or enhance 
the effects of aberration, depending on its orientation. For 
= — sin 6 cos aBpT /ri^c, the shapes of the phase portraits 
are summarised in Tableland in the Appendix. For a hollow 
cone, at those orientations where the pulse is double-peaked, 
the shapes are more complicated that those of a filled core. 
The tilt angles of the shapes are roughly proportional to 
the emission altitude and hence B^. The heart shapes and 
balloons seen in the Q-U plane also rotate about (U, Q) = 
(0,0) as the orientation and/or emission altitude changes. 

We close by summarising our main conclusions. 

(i) Stokes phase portraits contain additional information 
regarding a pulsar's orientation (q, i) and magnetic geome- 
try when used in conjunction with PA swings. 

(ii) The observed phase portraits of « 60% of the 24 
pulsars studied in this paper are asymmetric (Section 14. 3p . 
There are two possible causes for this: (1) the phase por- 
traits are distorted by relativistic aberration and are hence 
incompatible with a pure magnetic dipole at low emission 
altitudes, or (2) the axis of beam symmetry is tilted with re- 
spect to the magnetic axis, mimicing the effect of aberration 
even if the emission originates from a low altitude. 

(iii) If the former is true, the radio emission region in 
these objects is found to lie at r > O.lrLC, where the effects 
of relativistic aberration and a toroidal field are important. 
In some cases, one must have r > 0.2rLC, i-e. the emission 
comes from the outer magnetosphere . 

(iv) The idealised model of a pure or current-modified 
dipole at r > O.lrLC with a filled core beam and L — I cos 9o, 
or a hollow cone beam and L = /sin So, manages to ac- 



count for most gross features in the polarization data from 
all 26 pulsars studied in this paper, whether the pulse pro- 
files are single-peaked or double-peaked. Small longitudinal 
variations in / and L are needed to fit certain minor details 
(e.g. relative peak heights) in the observations. 

This paper is the first in a series. Future papers will 
examine the Stokes tomography of a force-free rotator, mil- 
lisecond pulsars, and the circularly polarized component of 
pulsar radio emission. 
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APPENDIX: GLOSSARY OF PHASE 
PORTRAIT DESCRIPTORS 

Figure [SS] gathers together, in a handy and easy-to-read ta- 
ble, the colloquial descriptors used to label phase portraits 
in the text. The parameters for each example are stated in 
the caption. 
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I sin Oq 


ESHII] 






Hollow cone 


I cos 0(j 
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Pure dipole 
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EHHH] 








I sin 9o 


ISHIS] 






Hollow cone 


I cos 9(1 
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Table 2. Summary of shapes seen in Stokes phase portraits for a pure dipole (Section |3)l. Examples of each shape appear in Figure [58l 
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Hollow cone trefoil, mosquito, twisted triangle {I-U) 
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balloon, heart, mosquito, interlocking ovals {Q-U) 
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Table 3. Summary of shapes seen in Stokes phase portraits for a current-modified dipole {Section^. Examples of each shape appear in 
Figure [58l 



L = I cos Oq 


L = I sin 9 Q 


Balloon, banana (I-Q) 


Balloon, banana {I-Q) 


Filled core Balloon, banana (I-U) 


Balloon, figure-eight (I-U) 


Balloon (Q-U) 


Balloon (Q-U) 


Balloon, mosquito, twisted triangle (I-Q) 
Hollow cone Balloon, banana, trefoil, mosquito (I-U) 
balloon, heart, figure-eight (Q-U) 


balloon, mosquito, figure-eight, trefoil (I-Q) 
Balloon, mosquito, trefoil (I-U) 
Balloon, heart, figure-eight (Q-U) 
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Figure 58. Colloquial descriptors of phase portraits encountered in this paper. Left column, top to bottom: (i) I-Q for pure dipole, filled 
core beam, (a,i) = (60°, 40°); (ii) U-Q for pure dipole, filled core beam, (ce,i) = (20°, 10°); (iii) I-U for pure dipole, filled core beam, 
(a,i) = (30°, 70°); (iv) I-U for pure dipole, hollow cone, (a,i) = (60°, 10°); (v) I-U for pure dipole, hollow cone, (a,i) = (50°, 30°). 
Right column, top to bottom: (i) I-U for current-modified dipole, hollow cone, (a,i) = (80°, 20°); (ii) I-U for pure dipole, hollow cone, 
(a,i) = (60°, 30°); (iii) I-U for pure dipole, hollow cone, (a,i) = (50°, 50°); (iv) Q-U for pure dipole, hollow cone, (a,i) = (70°, 40°); 
(v) I-Q for pure dipole, hollow cone, (a,i) = (50°, 40°). All examples share the same polarization pattern L/I = cosBq. 
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